Sign reversal of field-angle resolved heat capacity oscillations in a heavy fermion 
superconductor CeCoIn 5 and d x 2_ y 2 pairing symmetry 
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To identify the superconducting gap symmetry in CeCoIns (T c =2.3 K), we performed angle- 
resolved specific heat (C</>) measurements in a field rotated around the c-axis down to very low 
temperatures 0.05T C and detailed theoretical calculations. In a field of 1 T, a sign reversal of the 
fourfold angular oscillation in C<f, has been observed at T ~ 0.1T C on entering a quasiclassical 
regime where the maximum of corresponds to the antinodal direction, coinciding with the angle- 
resolved density of states (ADOS) calculation. The behavior, which exhibits minima along 
[110] directions, unambiguously allows us to conclude d x i_ y i symmetry of this system. The ADOS- 
quasiclassical region is confined to a narrow T and H domain within T/T c ~ 0.1 and 1.5 T (0.13if C 2)- 
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Investigation of superconductivity in strongly corre- 
lated electron systems has been a fascinating subject 
over the past few decades, since a realization of uncon- 
ventional non s-wave pairings can be expected due to 
strong electron-electron repulsion. A gap function of 
these unconventional superconductors (SCs) mostly has 
zeros (has nodes) along certain directions in the momen- 
tum space. The existence of nodes (point or line) can 
be inferred from power-law dependences in the physical 
quantities such as the specific heat (C) or the nuclear spin 
relaxation rate. To identify the gap structure is, however, 
a more formidable task. For instance, d xy and d x 2_ y 2 
symmetry differ from each other only in the position of 
line nodes (45° rotation of the latter becomes identical 
with the former). Phase-sensitive or direction-sensitive 
experiments are therefore needed to discriminate these 
gap structures. Very elegant phase-sensitive experiments 
have been done on the high-T c cuprates to establish pure 
d x 2_ y 2 symmetry, employing sophisticated tricrystal or 
tetracrystal geometries [1[ . 

Here we address the directional thermodynamic mea- 
surements. These techniques utilize that the low- lying 
quasiparticle excitations near the gap nodes become field 
orientation dependent [2|, |3[ in a low H and T region 
where the nodal quasiparticles contribute predominantly 
to thermodynamic quantities [1] ; The angle-resolved den- 
sity of states (ADOS) No((j)) at the Fermi energy becomes 
the largest (smallest) when H is parallel to the antinodal 
(nodal) directions. The nodal directions can be deter- 
mined from the angular variation of or the thermal 
conductivity in a rotating H (minima of or 
occur along the nodal directions). Thus those angle- 
resolved bulk thermodynamic experiments provide an in- 
dispensable spectroscopic method to locate the nodal di- 
rection. 

Up to now, the experiments have been performed 
on various anisotropic SCs to probe the gap struc- 



tures 0, @, 0, S, 0, Among these, we concen- 

trate here on the heavy fermion superconductor CeCoIns 
(T c =2.3 K) in which the gap structure has been con- 
troversial. The previous C<j> experiment (H rotated in 
the afr-plane) performed down to 0.3 K (T/T c >0.13) 
revealed a fourfold angular oscillation with the min- 
ima along [100], from which d xy symmetry was de- 
duced [6]. The same symmetry was also inferred from 
the anisotropy in the high- field superconducting phase 
of CeCoIns [ll[. This conclusion was, however, in dis- 
agreement with that of the k& measurements in which 
d x 2_ y 2 symmetry was derived [12]. Moreover, instability 
of the vortex lattice structure observed by small-angle 
neutron scattering experiments [H, Q has been shown 
to be consistent with the d x 2_ y 2 pairin g ITR [lo| . Point- 
contact Andreev reflection experiment [17[ and the s pin 
resonance observed in the superconducting state lH, 19[ 
also support d x 2_ y 2 symmetry in CeCoIns. 

Recently a resolution of this discrepancy has been put 
forward by Vekhter and co-workers [20|, [2l| . They the- 
oretically examined the phase diagram for of d-wave 
SCs in a wide H-T region and pointed out that apply- 
ing H along the gap nodes may result in maxima of 
in an inter mediate-T range. It is due to an inversion of 
the anisotropy in the finite energy DOS when the scat- 
tering of the quasiparticles off the vortices is considered. 
Thus, the ADOS-quasiclassical region in which C<j> takes 
minima for nodal direction may lie at still lower temper- 
atures. Such a crossover in with T has not yet been 
observed in nodal SCs. If this is the case for CeCoIn 5 , 
then the anisotropy of is expected to change the sign 
at a temperature c± 0.1T C [2(J,l2l|. In order to explore this 
possibility, we have extended the C<j> measurements on 
CeCoIns down to 0.05T C . Here we also perform more ac- 
curate calculations C<j> without resorting to the Pesch ap- 
proximation adopted by Vorontsov et al. [20[ and we also 
take into account Pauli paramagnetic effect and Fermi 
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FIG. 1: (a) C/T of CeCoIns at low temperatures measured 
in various fields applied along the [100] direction, (b) The 
electronic part (C — C n )/T obtained by subtracting the nu- 
clear contribution, (c) Field variation of the coefficient a(H) 
of the nuclear Schottky term. Dots and the solid line are the 
the experimental and the calculated results, respectively. 



FIG. 2: (color online). Left panel: Angular dependent specific 
heat data for selected fields and temperatures. Magnetic field 
is rotated in the afc-plane. Solid lines are fits to the expression 
Co+Cff(l — A4cos40). Right panel: Temperature variation 
of the normalized fourfold amplitude A4 obtained for H =3, 
2 and 1 T. 



velocity anisotropy effect, both of which are known to be 
important in this system [l3|, [TBj]. Those experimental 
and theoretical combined efforts pave a way to establish 
the angle-resolved thermodynamic measurements and 
^0 as important spectroscopic means to determine the 
pairing symmetry of unconventional superconductors. 

The single crystal of CeCoIn 5 (23 mg weight) used in 
the present experiment was a thin plate whose c axis 
oriented perpendicular to the largest face. measure- 
ments were done by a semiadiabatic heat pulse method. 
The sample was cooled using a dilution refrigerator (Ox- 
ford Kelvinox AST Minisorb) with a double sorption 
pump system built inside the insert (51 mm outer di- 
ameter) to allow continuous circulation of 3 He. Having 
no thick pumping tube outside, the insert could be eas- 
ily rotated using a stepper motor mounted at the top of 
a magnet Dewar. Horizontal fields were generated by a 
split-pair superconducting magnet. We confirmed that 
the addenda contribution, which was always less than 
10% of the sample specific heat, has no field-angle de- 
pendence. In a fixed field strength, we took the data 
every 2 degrees of the H rotation in the afr-plane cover- 
ing [100] through [010] directions. For each angle, the Cq 
value was determined by the average of three successive 
measurements. 

First we show in Fig. 1(a) the temperature variation 
of C/T obtained at H=0, 1, 2 and 3 T in a temper- 
ature range below 500 mK. At zero field, C/T above 



300 mK is linear to T as expected for the line node 
SC. At low temperatures below 200 mK, an upturn in 
C/T shows up. This upturn is due to a quadrupole 
splitting of the 115 In (1=9/2) and 59 Co (1=7/2) nuclear 
spins [H, . The nuclear contribution significantly in- 
creases by applying H and dominates the C/T data be- 
low 0.3 K. We subtracted the nuclear Schottky contribu- 
tion C n = a(H)/T 2 , where the coefficient a(H) is deter- 
mined at each field so that the resulting electronic contri- 
bution (C — C n )/T becomes linear to T at low tempera- 
tures as shown in Fig. 1(b) [24(. Dots in Fig. 1(c) denote 
a(H) obtained by this analysis, which show a H 2 like 
increase due to a nuclear Zeeman splitting. Of course, 
the above evaluation of a(H) might have some arbitrari- 
ness. Fortunately, the manner of subtracting C n does 
not strongly affect the analysis of Cq below, since C n 
is considered to be ^-independent as will be discussed 
later. The rapid increase of C n at low T, however, sets 
the lower bound of the temperature for the present Cq 
measurements to be ~0.12 K for H=l T and T >0.16 K 
for H=2 and 3 T, because of a small amplitude of the 
oscillation and a finite resolution (~0.5 %) of the specific 
heat measurements. 

In Fig. 2 (left panel), we show the nuclear subtracted 
C^ data for selected temperatures and fields. We al- 
ways observed a fourfold oscillation of Cq with no ap- 
preciable twofold part. We fit the data to an expression 
C — C n = Cq + Ctj(1 — A4 cos 40) , where Cq and Ch 
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are the zero field and field dependent parts of the elec- 
tronic specific heat, respectively, and <p is the angle of H 
measured from [100] direction. A4 is the normalized am- 
plitude of the fourfold oscillation which, for nodal SCs, 
becomes proportional to the zero-energy DOS anisotropy 
IA4I oc 1 - N (H || node) /N (H || antinode) in an ADOS- 
quasiclassical regime at low T and takes a nonzero value 
as H — > [3]. The A4 value thus obtained is shown in 
Fig. 2 (right panel) as a function of T for H=l, 2 and 
3 T. At T=300 mK, A4 is positive (C<f> minima along 
[100]) for all fields in agreement with the previous mea- 
surements [fjj]. For H=3 T, A4 is nearly T-independent 
below 300 mK. A4 for H=2 T shows a slight decrease on 
cooling, although it appears to remain positive as T — >0. 
When H is reduced to 1 T, A4 changes the sign around 
200 mK and becomes negative (C^ minima along [110]) 
at lower temperatures. 

In the evaluation of A4 above, we assumed that C n (H) 
is isotropic in the afr-plane. We confirmed the validity of 
this assumption by calculating C n (H) from the nuclear 
spin Hamiltonian [25] H n = (/w/q/6)[3/| - 1(1 + 1) + 
rj(I% — Iy)] + 7^/(1 + K)H, where vq and n are the 
parameters describing a quadrupole splitting, and K is 
the Knight shift tensor. In CeCoIns, there are four low- 
symmetry (77 7^ 0) In sites (per unit cell) on the lateral 
faces of the unit cell. K also has an in-plane anisotropy in 
these In sites. Each site of these then produces a twofold 
in-plane anisotropy in C n (H). When averaged over the 
four sites, however, the twofold anisotropy is cancelled, 
leaving no (fourfold) <\> dependence in C n (H). The solid 
line in Fig. 1(c) is the field variation of a(H) calculated 
with the parameters taken from Ref. |25[ , which is in good 
agreement with the experimentally evaluated results. We 
may therefore conclude that the behavior of A±(T,H) 
above is the intrinsic properties in the superconducting 
state. 

In Fig. 3, we show the contour plot of A*(T, if) ob- 
tained by the present measurements. Our present results 
are fully compatible with the previous experiment 
performed in the region T > 0.1 3T C and H > 0.17i^ C 2, 
where A±(T,H) is positive and takes minima along 
[100] directions. The sign reversal in A±(T,H) is ob- 
served only at much lower temperatures (T < 0.1 2T C ) 
and lower fields (H < 0.15H C 2). As will be shown in the 
following, the crossover temperature of the sign rever- 
sal (~0.12T C ) well agrees with the theoretical prediction. 
These observations indicate that, for the first time, the 
low temperature ADOS-quasiclassical region is reached 
experimentally. In this region, takes minima along 
[110] directions implying unambiguously that the gap 
symmetry of CeCoIns is d x 2_ y 2. 

In order to understand the above data of C<j> or A4 
and also to help establishing the angle-resolved thermo- 
dynamic measurements, we calculate the angle-resolved 
specific heat C^T^H) accurately in the quasiclassical 
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FIG. 3: Contour plot of the normalized fourfold amplitude 
A±(T,H) of CeCoIns obtained by the present experiment. T 
and H are normalized by T c =2.3 K and if C 2=H-5 T, respec- 
tively. 

Eilenberger theory without resorting to the Pesch ap- 
proximation adopted by Vorontsov et al [20[ . The Eilen- 
berger equations are read as 

{uj n + inB + v F • [V + iA(r)]}f = A(r)4>(k)g (1) 
{u n + ipB - v F • [V - iA(r)]}p = A*(r)4>(k)g (2) 

where the quasiclassical Green functions are related to 
g = (1 — fP) 1 / 2 . The d x 2_ y 2 pairing function is 
given by 0(fc) = cos 2$ where is azimuthal angle 
in /c-space. The Pauli paramagnetic parameter \i con- 
trols its strength. In the rippled cylindrical Fermi sur- 
face, oc (v cos v sin v z sin k z ) with v z = 0.5 and 
v = 1 + /3 cos AO where /3 signifies the Fermi velocity 
anisotropy. The anisotropy ratio of the coherent lengths 
ic/iab ~ 0.5 by this v z . The notation and the proce- 
dure of the calculation are the same as before [26[ al- 
beit the field direction here is in the afr-plane. After 
solving those equations self-consistently, we evaluate the 
Helmholtz free energy and the entropy, which lead to 
the specific heat and the fourfold oscillation amplitude 
MT,H). 

We have done extensive numerical computations of 
A±(T,H) for various \i and /3 values, including the pre- 
vious case with \i = and (3 = whose result is basi- 
cally consistent with that in Ref. [2(j. Here we take into 
account the Pauli paramagnetic effect and the Fermi ve- 
locity anisotropy effect, both of which are known to be 
important in CeCoIns. In Fig. 4 we show the contour 
plot of A A (T,H) for /i = 2 and (3 = 0.5. It is seen 
that (i) The ADOS-quasiclassical region is confined to 
a narrow region in low H and T bounded by H s and 
T s . This region is characterized by A4 < where the 
oscillation sense coincides with the angle-resolved DOS, 
namely the maximum occurs when the field direction is 
along the antinodal [100] direction expected by d x 2_ y 2 
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ADOS-quasiclassical region where the oscillation maxi- 
mum coincides with the anti- nodal direction [100]. Ex- 
cept for trivial sign changes near H C 2 (see for example 
@), this is the first time. Thus the gap symmetry in 
CeCoIns is d x 2_ y 2. Those findings also help establish- 
ing the angle-resolved thermodynamic measurements as 
a unique and indispensable spectroscopic method to iden- 
tify the gap location, or the pairing symmetry. 

This work has partly been supported by a Grant-in- 
Aid for Scientific Research on Innovative Areas "Heavy 
Electrons" (No. 20102007) of the Ministry of Education, 
Culture, Sports and Technology, Japan. 



FIG. 4: Topography of the specific heat oscillation amplitude 
A4 (in units of percent) calculated under the paramagnetic 
effect parameter \i — 2.0, the Fermi velocity anisotropy (3 = 
0.5, and GL parameter k = 89. The inside of the curved area 
around the origin is the ADOS-quasiclassical region bounded 
by H s and T s . The calculations are done for 6 x 7 grid points 
in (T, H) and interpolated. 



symmetry, (ii) The other region with A4 > in Fig. 
4 is occupied by the reversed oscillation region. The 
maximum oscillation is located around T/T c ~ 0.25 and 
H/H c2 ~ 0.35. (hi) The overall landscape A 4 (T,H) well 
explains the experimental data in Fig. 3. Thus we firmly 
conclude that the newly discovered sign reversal is un- 
derstood as d x 2_ y 2 symmetry, (iv) We note the fallow- 
ings: The ADOS-quasiclassical region in Fig. 4 ends at 
H S /H C 2 ~ 0.3 and T s /T c ~ 0.1 while the experimental 
data show H S /H C 2 ~ 0.15 and T s /T c ~ 0.1. There is a 
large discrepancy in H s value. Note that H S /H C 2 ~ 0.4 
for the d-wave case with /i = (3 = and H S /H C 2 ~ 0.35 
for the d-wave case with n = 2 and (3 = 0. Thus the up- 
per bound H s is sensitive for material parameters while 
T s /T c is rather independent of those. This robustness 
is important when interpreting data, namely to find the 
sign change the T sweep is more effective than H sweep. 

Thermal conductivity experiment by Izawa et 
al. 12] on CeCoIn 5 is done outside the ADOS- 



quasiclassical region (T/T c > 0.2). They correctly con- 
clude the d x 2_ y 2 symmetry, which is fully consistent with 
the present identification. We comment on the recent 
experiment by Park et al. [To[ on CeRhlns under pressure 
taken at H/H c2 ~ 0.13 and T/T c = 0.14. The data show 
the maximum at [110] direction. Since that temperature 
is just at the boundary of the ADOS-quasiclassical re- 
gion, we should be careful to conclude the gap symmetry. 
Based on the present result, we suggest an experiment by 
changing T under fixed field to check whether the data 
is inside or outside that region. 

In summary, we have found the sign change of the spe- 
cific heat oscillation and discovered the (i^, T) region, or 
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